Genetic-Algorithm-Based Inverse Optimization Identification Method for Hot-Temperature Constitutive Model Parameters of Ti6Al4V Alloy

A precise constitutive model is the foundation and key to finite element simulation in material volume forming and the optimization of the hot working process. Hence, to build a precise constitutive model, a method based on a genetic algorithm (GA) for the inverse optimization identification of parameters is presented in this paper. The idea of this method is to continuously adjust the model parameters through GA until the objective function reaches the minimum value. In this study, hot compression experiments were performed on the Gleeble-1500D thermal simulator at temperatures ranging from 800 °C to 1000 °C and strain rates of 0.01 s−1 to 1 s−1. The Arrhenius-type (A-T) model considering strain compensation and the Johnson–Cook (JC) model considering the coupling effects of strain, temperature and strain rate were constructed, respectively, by using the regression method and the parameter inverse optimization identification method. For the purposes of comparing and verifying the reliability of the predictions of the two established constitutive models, the correlation coefficient (R), average absolute relative error (AARE), and relative error (RE) were adopted. The results show that both the optimized A-T model and the optimized JC model have high prediction accuracy. Compared to the optimized JC model, the optimized A-T model demonstrated a higher correlation coefficient, by 0.003, and a lower average absolute relative error, by 1.43%. Furthermore, the relative error distribution of the optimized A-T model was found to be more concentrated than that of the optimized JC model. These results suggest that the A-T model is more appropriate than the JC model for characterizing the high-temperature deformation behavior of Ti6Al4V alloy.


Introduction
Because of the outstanding characteristics of high specific strength, corrosion resistance, and biocompatibility, titanium alloys have found extensive applications in various fields such as aerospace, the chemical industry, automobiles, and medical implants [1][2][3]. Among all titanium alloys, Ti6Al4V alloy is the most frequently used high-strength duplex (α+β) titanium alloy [4]. At room temperature, the Ti6Al4V alloy has poor formability due to its densely packed hexagonal crystal structure, which has only three slip systems. This makes it difficult to produce complex-shaped parts [5]. Thus, the hot formation of Ti6Al4V alloy has always been the focus of research by scholars [6,7]. Different hot deformation conditions will cause different microstructure evolutions of the Ti6Al4V alloy during hot forming, resulting in a highly nonlinear relationship between stress and temperature, strain rate, and strain at the macroscopic level [8]. Therefore, studying the hot deformation behavior of metals and establishing constitutive models that precisely characterize the ε ·(1 + ε) m 5 T ·e m 7 ε · . ε m 8 T ·T m 9 (1) For the solution of parameters m 3 and m 8 , the temperature and strain are firstly controlled to be constant. Then, the terms of the above formula for temperature and strain are all constants, which are set as K. Equation (1) can be changed to: ε m 3 +m 8 T ·K (2) Equation (2) can be transformed by taking the natural logarithm of both sides: From Equation (3) we can obtain: The average slope of the fitted lines based on the relationship between lnσ and ln . ε can be used to determine the value of m 3 + m 8 T. Then, the relationship between m 3 + m 8 T and T can be linearly fitted. It can be clearly seen that the value of m 8 is given by the slope of the fitted line, while the value of m 3 is given by the intercept. Likewise, the remaining parameters of the Hansel-Spittel constitutive model can also be derived using this approach, if the maximum and minimum values of the slope of the fitted straight lines are quite different. At this time, the selected average value will cause a large error in m 3 + m 8 T. The same goes for other parameters. Finally, the accuracy of the model cannot be guaranteed. The Hansel-Spittel model established by Chadha [31] had a maximum error of 14% between the predicted values and experimental values, and at low strain rates the model's ability to predict the alloy's softening behavior was limited. Consequently, the constitutive models established through the regression method often fail to precisely predict the flow behavior of alloys.
In recent years, the inverse analysis method has been widely used in the identification process of parameters in various industries [33][34][35][36]. For example, Levasseur S [33] used the inverse analysis method to identify the parameters of the soil constitutive model. Gajewski T [34] quickly calibrated a complete set of parameters for concrete using an inverse analysis method. Among them, the genetic algorithm (GA) is a widely used global optimization algorithm. It seeks the optimal solution to a problem by imitating the phenomena of replication, crossover, and mutation in natural selection and heredity. The main method of using GA to solve the optimal solution of parameters is through iterative operation. Some optimization algorithms are prone to become trapped in local minima and can exhibit the "dead loop" phenomenon, but GA can overcome this shortcoming very well. Compared with traditional optimization methods (enumeration, heuristic, etc.), GA has less calculation time and characteristics of good convergence and robustness. Both Wu [37] and Chen W [38] employed the genetic algorithm (GA) to perform inverse optimization of the material parameters in the A-T model. It was found that the A-T model's forecasting accuracy greatly increased. Through investigation, it was found that there are few studies on the optimization of constitutive model parameters for the Ti6Al4V alloy. Therefore, a parameter inverse optimization identification method based on GA to develop the constitutive model of Ti6Al4V alloy was chosen in this paper.
This paper aims to develop a precise constitutive model and lay the foundation for the subsequent reasonable design of Ti6Al4V thermal processing parameters and finite element simulation. The Ti6Al4V alloy's hot deformation characteristics were analyzed using the Gleeble-1500D thermal simulator. In order to accurately identify the parameters of the constitutive model, a parameter inverse optimization identification method based on the genetic algorithm was proposed in this paper. The objective function was the accumulated error between the stress calculated by the constitutive model and the experimental stress. This method iteratively adjusts the parameters of the constitutive model through the genetic algorithm until the objective function reaches the minimum value, and the output parameters are the optimal solution of the parameters of the constitutive model. Finally, the prediction accuracy of the two models was quantitatively evaluated using standard statistical parameters, and the constitutive model with the highest prediction accuracy was selected.

Experimental Method
The chemical composition of the Ti6Al4V alloy can be seen in Table 1. To investigate the hot deformation characteristics of the Ti6Al4V alloy and establish its constitutive model, thermal compression tests were carried out on the Gleeble-1500D thermal simulator. The size of the compression samples was φ8 mm × 12 mm. To make the deformation of the samples uniform and mitigate or eliminate the impact of the frictional forces on the experimental results, graphite was uniformly applied to both ends of the cylinder in contact with the indenter of the testing machine. The deformation temperatures were 800 • C, 850 • C, 900 • C, 950 • C, and 1000 • C; the strain rates were 0.01 s −1 , 0.1 s −1 , and 1 s −1 . The thermal compression experiment process is shown in Figure 1. Firstly, a heating rate of 10 • C per second was applied to each cylindrical sample. After reaching the predetermined temperature, it was kept warm for 180 s to remove the internal temperature gradient of the samples. Lastly, compression tests were conducted and each sample was compressed by 50%. This paper aims to develop a precise constitutive model and lay the foundation for the subsequent reasonable design of Ti6Al4V thermal processing parameters and finite element simulation. The Ti6Al4V alloy's hot deformation characteristics were analyzed using the Gleeble-1500D thermal simulator. In order to accurately identify the parameters of the constitutive model, a parameter inverse optimization identification method based on the genetic algorithm was proposed in this paper. The objective function was the accumulated error between the stress calculated by the constitutive model and the experimental stress. This method iteratively adjusts the parameters of the constitutive model through the genetic algorithm until the objective function reaches the minimum value, and the output parameters are the optimal solution of the parameters of the constitutive model. Finally, the prediction accuracy of the two models was quantitatively evaluated using standard statistical parameters, and the constitutive model with the highest prediction accuracy was selected.

Experimental Method
The chemical composition of the Ti6Al4V alloy can be seen in Table 1. To investigate the hot deformation characteristics of the Ti6Al4V alloy and establish its constitutive model, thermal compression tests were carried out on the Gleeble-1500D thermal simulator. The size of the compression samples was ϕ8 mm × 12 mm. To make the deformation of the samples uniform and mitigate or eliminate the impact of the frictional forces on the experimental results, graphite was uniformly applied to both ends of the cylinder in contact with the indenter of the testing machine. The deformation temperatures were 800 °C, 850 °C, 900 °C, 950 °C, and 1000 °C; the strain rates were 0.01 s −1 , 0.1 s −1 , and 1 s −1 . The thermal compression experiment process is shown in Figure 1. Firstly, a heating rate of 10 °C per second was applied to each cylindrical sample. After reaching the predetermined temperature, it was kept warm for 180 s to remove the internal temperature gradient of the samples. Lastly, compression tests were conducted and each sample was compressed by 50%.

Analysis of Experimental Results of the Ti6Al4V Alloy
The true stress-strain curves of the Ti6Al4V alloy under different deformation conditions can be observed in Figure 2. Holding the strain rate constant, an increase in temperature resulted in a decrease in flow stress. From Figure 2c, it can be observed that the material's peak stress decreased from 336.95 MPa at 800 °C to 71.09 MPa at 1000 °C. This

Analysis of Experimental Results of the Ti6Al4V Alloy
The true stress-strain curves of the Ti6Al4V alloy under different deformation conditions can be observed in Figure 2. Holding the strain rate constant, an increase in temperature resulted in a decrease in flow stress. From Figure 2c, it can be observed that the material's peak stress decreased from 336.95 MPa at 800 • C to 71.09 MPa at 1000 • C. This was because the thermal activation energy of the alloy increased with an increase in deformation temperature, leading to a more pronounced softening effect. At a constant temperature, the flow stress rose with an increase in the strain rate. According to Figure 2d, the peak stress of the material increased from 72.26 MPa at 0.01 s −1 to 165.47 MPa at 1 s −1 .This is because the strain rate was too fast, and dislocations accumulated rapidly, resulting in the obvious work hardening effect of the Ti6Al4V alloy.
Ti6Al4V alloy was α phase. The α phase is a hexagonal close-packed structure wi stacking fault energy, which makes it relatively easy for dislocations to move in th tal, and dynamic recrystallization occurs easily. Therefore, the dynamic recrystalli softening phenomenon could be clearly observed at 800-900 °C. At the 800 °C/0.1 formation condition, the flow stress decreased by 82.98 MPa from its peak value, as in Figure 2b. The dynamic recovery phenomenon appeared under the deformation tion of 900 °C/1 s −1 , as shown in Figure 2c. This is because the recrystallization tim relatively short at high strain rates, so the stress rose slowly to 165.47 MPa after q reaching 152.03 MPa. The α phase will gradually transform into the β phase above 850 °C. In cont the α phase, the β phase is a body-centered cubic structure with high stacking fault e and many slip systems, resulting in insufficient deformation energy to provide the d force for dynamic recrystallization. Therefore, above 950 °C, dynamic recovery is th softening mechanism of the alloy, and the stress has no obvious downward trend. In addition, since the temperature was below 850 • C, the main microstructure of the Ti6Al4V alloy was α phase. The α phase is a hexagonal close-packed structure with low stacking fault energy, which makes it relatively easy for dislocations to move in the crystal, and dynamic recrystallization occurs easily. Therefore, the dynamic recrystallization softening phenomenon could be clearly observed at 800-900 • C. At the 800 • C/0.1 s −1 deformation condition, the flow stress decreased by 82.98 MPa from its peak value, as shown in Figure 2b. The dynamic recovery phenomenon appeared under the deformation condition of 900 • C/1 s −1 , as shown in Figure 2c. This is because the recrystallization time was relatively short at high strain rates, so the stress rose slowly to 165.47 MPa after quickly reaching 152.03 MPa.
The α phase will gradually transform into the β phase above 850 • C. In contrast to the α phase, the β phase is a body-centered cubic structure with high stacking fault energy and many slip systems, resulting in insufficient deformation energy to provide the driving force for dynamic recrystallization. Therefore, above 950 • C, dynamic recovery is the main softening mechanism of the alloy, and the stress has no obvious downward trend.

Parameter Inverse Optimization Identification Method Based on GA
The idea of the parameter inverse optimization identification method is to continuously adjust the parameters of the constitutive model through GA. When the objective function approaches zero infinitely, the output parameters are the optimal solution. Figure 3 shows the logical flowchart of this method. Firstly, input the data obtained from thermal compression experiments and the expression for the constitutive model. After that, the GA is invoked. Finally, set the initial value and iteration range of the constitutive model parameters.

Parameter Inverse Optimization Identification Method Based on GA
The idea of the parameter inverse optimization identification method is to continuously adjust the parameters of the constitutive model through GA. When the objective function approaches zero infinitely, the output parameters are the optimal solution. Figure 3 shows the logical flowchart of this method. Firstly, input the data obtained from thermal compression experiments and the expression for the constitutive model. After that, the GA is invoked. Finally, set the initial value and iteration range of the constitutive model parameters.  The material parameters obtained by the regression method were set as the initial values of the GA iteration. Sections 3.2 and 3.3 describe how to determine constitutive model parameters using the regression method. The setting of the model parameters' iteration range refers to the initial values. In this paper, the cumulative error between the stress predicted by the Ti6Al4V alloy constitutive model and the experimentally measured stress is set as the objective function, as shown in Formula (5): where exp i σ is the i-th stress measured experimentally and cal i σ is the stress of the i-th data calculated by the constitutive model.

Establishment of A-T Constitutive Model for the Ti6Al4V Alloy
An exponential equation that includes the Zener-Hollomon ( Z ) parameter is commonly used to describe the effect of temperature and strain rate on the thermal deformation behavior of metal materials [39]: Furthermore, the Z parameter can be described by various functions in distinct stress ranges [40]. The expressions for each range are provided below: The material parameters obtained by the regression method were set as the initial values of the GA iteration. Sections 3.2 and 3.3 describe how to determine constitutive model parameters using the regression method. The setting of the model parameters' iteration range refers to the initial values. In this paper, the cumulative error between the stress predicted by the Ti6Al4V alloy constitutive model and the experimentally measured stress is set as the objective function, as shown in Formula (5): where σ i exp is the i-th stress measured experimentally and σ i cal is the stress of the i-th data calculated by the constitutive model.

Establishment of A-T Constitutive Model for the Ti6Al4V Alloy
An exponential equation that includes the Zener-Hollomon (Z) parameter is commonly used to describe the effect of temperature and strain rate on the thermal deformation behavior of metal materials [39]: Furthermore, the Z parameter can be described by various functions in distinct stress ranges [40]. The expressions for each range are provided below: The characters in Formulas (6)- (8) are explained in Table 2. The solution formula of flow stress can be obtained by combining Equations (6) and (7) and the third equation in Equation (8), and the expression is as follows:

Parameters Description
True stress (MPa) A, n, α, β, and n Material parameters Furthermore, combining Equations (6)-(8), the following relationship can be derived: The natural logarithm of both sides of the three equations in Formula (10) can be taken to yield: In the equations, B and C are material parameters. The impact of strain was factored in when determining the values of material parameters. In this paper, the relationship between strain and material parameters was constructed using polynomial regression. The material parameters were determined using experimental data within the true strain range of 0.05-0.65. The sampling interval of strain was 0.05. Here, the process for determining the material parameters is illustrated using the example of a strain of 0.1.
When the strain is 0.1, the experimental data are substituted into the first two equations of Equation (11). The material parameters β and n can be determined using a linear regression approach. In Figure 4a,b, were parameters are determined by computing the inverse of the slope of the fitted lines. The data analysis revealed that the two parameters had average values of 0.062997 MPa −1 and 5.938013, respectively. The material parameter α could be derived as 0.010609 MPa −1 using the equation α = β/n .
Similarly, using the results of fitting straight lines in Figure 5a,b, the parameters n, Q, and ln A were calculated to have values of 3.839814, 630,040.1 J·mol −1 , and 60.44771 s −1 , respectively.
Likewise, parameters under other strains were determined in the same way. The results can be found in Table 3.  Likewise, parameters under other strains were determined in the same way. The results can be found in Table 3.    Likewise, parameters under other strains were determined in the same way. The results can be found in Table 3.  Using the parameter values at 0.1 strain in Table 3 as the initial values for the GA iteration and executing the flowchart in Section 3.1, under this strain, the change of the objective function value of the A-T model with the number of iterations is shown in Figure 6. As the iterations proceeded, the objective function value eventually converged to around   Table 3 as the initial values for t iteration and executing the flowchart in Section 3.1, under this strain, the change objective function value of the A-T model with the number of iterations is shown in 6. As the iterations proceeded, the objective function value eventually conver around 0.00289. The output of the logic diagram was the optimal solution of the co tive model parameters under the strain. Similarly, using the same method, the o solutions of material parameters under other strains are shown in Table 4.   To improve the accuracy of polynomial regression, sixth-order polynomial regression was performed on the material parameters and true strain: The relationship between true strain and material parameters is depicted in Figure 7. The four parameters changed greatly after optimization. tables 5 and 6 are the sixthorder polynomial coefficient values of the unoptimized A-T model and the optimized A-T model, respectively. ln ( ) The relationship between true strain and material parameters is depicted in Figure 7. The four parameters changed greatly after optimization. Tables 5 and 6 are the sixth-order polynomial coefficient values of the unoptimized A-T model and the optimized A-T model, respectively.

Establishment of the JC Constitutive Model for the Ti6Al4V Alloy
Apart from the A-T model, the JC model is also extensively employed during metal forming processes owing to its few parameters and simple fitting process [41]. In this paper, a modified JC constitutive model is proposed. Its expression is: The meanings expressed by each character in Formula (13) are listed in Table 7. In this paper, . ε re f is equivalent to 0.01 s −1 and T re f is denoted as 800 • C. The melting point of the Ti6Al4V alloy is 1660 • C [42]. Table 7. Description of parameters in the JC model.

Polynomial coefficients
The flowchart for determining the material parameters of the JC model is presented in Figure 8. Firstly, under the condition of 800 • C/0.01 s −1 , the relationship between ε and σ is fitted by a polynomial. The polynomial coefficients correspond to the values of parameters A 0 , A 1 , A 2 , and A 3 , consecutively.
Then, the relationship between σ/A(ε) and ln . ε * was fitted linearly at 800 • C, and the values of B(ε) at different strains could be determined by the slopes of the fitted lines. Then, fit the relationship between B(ε) and ε through the third-order polynomial. Similar to the previous scenario, the coefficients of the third-order polynomial are indicative of the values of parameters B 0 , B 1 , B 2 , and B 3 , in that order. Finally, a linear fit was made to the relationship between ln and T * . The values of C(ε) + D(ε) ln . ε * corresponded to the slopes of the fitted lines. When the strain rate was fixed at 0.1 s −1 , determining the value of parameters C(ε) at various strains was a simple task, and calculating parameters D(ε) was also possible. Likewise, by employing identical third-order polynomial regression, the parameters C 0 − C 3 and D 0 − D 3 could be determined. In summary, the JC model's parameters are listed in Table 8. These parameters were employed as the initial values for the GA iteration. Figure 9 is the optimization curve of the objective function of the JC model. The value of the objective function gradually decreased with the number of iterations, and finally tended to be around 0.00373. At this stage, the output results represent the optimal solution for the JC model parameters, as listed in Table 9.

2, 3) Polynomial coefficients
The flowchart for determining the material parameters of the JC model is prese in Figure 8. Firstly, under the condition of 800 °C/0.01 s −1 , the relationship betwee and  is fitted by a polynomial. The polynomial coefficients correspond to the valu     When the strain rate was fixed at 0.1 s −1 , determining the value of parameters various strains was a simple task, and calculating parameters D(ε) was also p Likewise, by employing identical third-order polynomial regression, the par  Table 8. These paramet employed as the initial values for the GA iteration. Figure 9 is the optimization the objective function of the JC model. The value of the objective function gradu creased with the number of iterations, and finally tended to be around 0.00373 stage, the output results represent the optimal solution for the JC model param listed in Table 9. Table 8. The values of third-order polynomial coefficients for the unoptimized JC model.

Verification and Comparison of Prediction Accuracy of Two Constitutive Models for the Ti6Al4V Alloy
The predicted values of the constitutive models established in the previous section were compared to the experimental values. The results are depicted in Figure 10. For the A-T model, the optimization effect was evident at 800 • C. For the JC model, the optimization effect was more pronounced under the three deformation conditions of 900 • C/1 s −1 , 950 • C/1 s −1 , and 1000 • C/1 s −1 , as depicted in Figure 10f. Furthermore, since 800 • C/0.01 s −1 was the reference condition of the JC model, therefore the JC model exhibited higher predictive reliability than the A-T model under this deformation condition.

Verification and Comparison of Prediction Accuracy of Two Constitutive Models for the Ti6Al4V Alloy
The predicted values of the constitutive models established in the previous section were compared to the experimental values. The results are depicted in Figure 10. For the A-T model, the optimization effect was evident at 800 °C. For the JC model, the optimization effect was more pronounced under the three deformation conditions of 900 °C/1 s −1 , 950 °C/1 s −1 , and 1000 °C/1 s −1 , as depicted in Figure 10f. Furthermore, since 800 °C/0.01 s −1 was the reference condition of the JC model, therefore the JC model exhibited higher predictive reliability than the A-T model under this deformation condition. To quantitatively assess the reliability of the above two models, R and AARE were introduced. Their formulas can be expressed as: To quantitatively assess the reliability of the above two models, R and AARE were introduced. Their formulas can be expressed as: In the formula, σ e is the experimentally measured stress, σ cal is the stress calculated by the model, σ e is the average of all σ e measured in the experiment, and σ cal is the average of all σ cal .
The value of R indicates the degree of correlation between the predicted values and the experimental values. However, because R does not consider bias, higher R values do not always mean a better predictiveness of the constitutive model. AARE is an unbiased statistical parameter used to measure the predictiveness of a model or equation. The smaller the AARE value, the higher the predictiveness of the model.
It can be seen from Figure 11 that the optimized A-T model has an R value of 0.9977 and an AARE value of 4.92%. The R value of the optimized A-T model is 0.0086 greater than that of the unoptimized A-T model, and the AARE value of the optimized A-T model is 2.98% lower than that of the unoptimized A-T model. The optimized JC model has an R value of 0.9947 and an AARE value of 6.35%. Similarly, the optimized JC model showed an increase in R value and a decrease in AARE value, with an increase of 0.0172 and a decrease of 4.81%, respectively. This shows that the proposed parameter inverse optimization identification method is feasible.  The value of R indicates the degree of correlation between the predicted values and the experimental values. However, because R does not consider bias, higher R values do not always mean a better predictiveness of the constitutive model. AARE is an unbiased statistical parameter used to measure the predictiveness of a model or equation. The smaller the AARE value, the higher the predictiveness of the model.
It can be seen from Figure 11 that the optimized A-T model has an R value of 0.9977 and an AARE value of 4.92%. The R value of the optimized A-T model is 0.0086 greater than that of the unoptimized A-T model, and the AARE value of the optimized A-T model is 2.98% lower than that of the unoptimized A-T model. The optimized JC model has an R value of 0.9947 and an AARE value of 6.35%. Similarly, the optimized JC model showed an increase in R value and a decrease in AARE value, with an increase of 0.0172 and a decrease of 4.81%, respectively. This shows that the proposed parameter inverse optimization identification method is feasible.
As shown in Figure 11a,c, after optimization, both the A-T model and the JC model had high predictability for the Ti6Al4V alloy. After comparing the A-T model and the JC model, it was found that the R value of the former was 0.003 higher than that of the latter, while the AARE value of the former was 1.43% lower than that of the latter. It was shown that, after optimization, the A-T model was better suited than the JC model for describing the high-temperature deformation behavior of the Ti6Al4V Alloy.  Figure 11a,c, after optimization, both the A-T model and the JC model had high predictability for the Ti6Al4V alloy. After comparing the A-T model and the JC model, it was found that the R value of the former was 0.003 higher than that of the latter, while the AARE value of the former was 1.43% lower than that of the latter. It was shown that, after optimization, the A-T model was better suited than the JC model for describing the high-temperature deformation behavior of the Ti6Al4V Alloy.

As shown in
For the purpose of further assessing the prediction accuracy of the two optimized models, the RE between experimental stress and model-calculated stress was analyzed. The formula for RE is as follows: In the formula, the definitions of σ e and σ cal are the same as in the previous formulas. Figure 12 shows the RE distributions of the optimized A-T model and the optimized JC model. It can be seen from Figure 12 that the RE fluctuation range of the optimized A-T model was between −20% and 20%, while the maximum RE of the optimized JC model reached −35%. The mathematical expectation and standard deviation of the optimized A-T model RE were 0.67406 and 6.34041, respectively, and those of the optimized JC model were −2.36206 and 9.77954, respectively, as determined by calculation. This indicates that the RE distribution of the optimized A-T model was more concentrated and tended to a smaller RE value. This further confirms that the A-T model had a greater predictive accuracy compared to the JC model. For the purpose of further assessing the prediction accuracy of the two optimized models, the RE between experimental stress and model-calculated stress was analyzed. The formula for RE is as follows: In the formula, the definitions of e  and  cal are the same as in the previous formulas. Figure 12 shows the RE distributions of the optimized A-T model and the optimized JC model. It can be seen from Figure 12 that the RE fluctuation range of the optimized A-T model was between −20% and 20%, while the maximum RE of the optimized JC model reached −35%. The mathematical expectation and standard deviation of the optimized A-T model RE were 0.67406 and 6.34041, respectively, and those of the optimized JC model were −2.36206 and 9.77954, respectively, as determined by calculation. This indicates that the RE distribution of the optimized A-T model was more concentrated and tended to a smaller RE value. This further confirms that the A-T model had a greater predictive accuracy compared to the JC model.

Conclusions
In this paper, the thermal deformation behavior of Ti6Al4V alloy was analyzed. The A-T model and the JC model were, respectively, established by using the parameter inverse optimization identification method. Finally, the accuracy of the two models' predictions was quantitatively compared by introducing standard statistical parameters. The specific conclusions obtained are: (1) The true stress of Ti6Al4V alloy increases with the increase of strain rate and decreases with the increase of temperature. Below 850 °C, the Ti6Al4V alloy will show obvious dynamic recrystallization characteristics. Since the α phase will gradually transform into the β phase above 850 °C, and the β phase has high stacking fault energy, the flow stress shows a dynamic recovery phenomenon above 950 °C. (2) The A-T model and JC model established by the inverse optimization identification method both exhibit higher R values and smaller AARE values, demonstrating the feasibility of the method proposed in this paper.

Conclusions
In this paper, the thermal deformation behavior of Ti6Al4V alloy was analyzed. The A-T model and the JC model were, respectively, established by using the parameter inverse optimization identification method. Finally, the accuracy of the two models' predictions was quantitatively compared by introducing standard statistical parameters. The specific conclusions obtained are: (1) The true stress of Ti6Al4V alloy increases with the increase of strain rate and decreases with the increase of temperature. Below 850 • C, the Ti6Al4V alloy will show obvious dynamic recrystallization characteristics. Since the α phase will gradually transform into the β phase above 850 • C, and the β phase has high stacking fault energy, the flow stress shows a dynamic recovery phenomenon above 950 • C.

Institutional Review Board Statement:
The study did not require ethical approval. Studies are not involving humans or animals.
Informed Consent Statement: "Not applicable" for studies not involving humans.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to these data being part of ongoing research.

Conflicts of Interest:
The authors declare no conflict of interest.